RNAstructure: software for RNA secondary structure prediction and analysis

<p>Abstract</p> <p>Background</p> <p>To understand an RNA sequence's mechanism of action, the structure must be known. Furthermore, target RNA structure is an important consideration in the design of small interfering RNAs and antisense DNA oligonucleotides. RNA secondary structure prediction, using thermodynamics, can be used to develop hypotheses about the structure of an RNA sequence.</p> <p>Results</p> <p>RNAstructure is a software package for RNA secondary structure prediction and analysis. It uses thermodynamics and utilizes the most recent set of nearest neighbor parameters from the Turner group. It includes methods for secondary structure prediction (using several algorithms), prediction of base pair probabilities, bimolecular structure prediction, and prediction of a structure common to two sequences. This contribution describes new extensions to the package, including a library of C++ classes for incorporation into other programs, a user-friendly graphical user interface written in JAVA, and new Unix-style text interfaces. The original graphical user interface for Microsoft Windows is still maintained.</p> <p>Conclusion</p> <p>The extensions to RNAstructure serve to make RNA secondary structure prediction user-friendly. The package is available for download from the Mathews lab homepage at <url>http://rna.urmc.rochester.edu/RNAstructure.html</url>.</p

From Structure Prediction to Genomic Screens for Novel Non-Coding RNAs

Non-coding RNAs (ncRNAs) are receiving more and more attention not only as an abundant class of genes, but also as regulatory structural elements (some located in mRNAs). A key feature of RNA function is its structure. Computational methods were developed early for folding and prediction of RNA structure with the aim of assisting in functional analysis. With the discovery of more and more ncRNAs, it has become clear that a large fraction of these are highly structured. Interestingly, a large part of the structure is comprised of regular Watson-Crick and GU wobble base pairs. This and the increased amount of available genomes have made it possible to employ structure-based methods for genomic screens. The field has moved from folding prediction of single sequences to computational screens for ncRNAs in genomic sequence using the RNA structure as the main characteristic feature. Whereas early methods focused on energy-directed folding of single sequences, comparative analysis based on structure preserving changes of base pairs has been efficient in improving accuracy, and today this constitutes a key component in genomic screens. Here, we cover the basic principles of RNA folding and touch upon some of the concepts in current methods that have been applied in genomic screens for de novo RNA structures in searches for novel ncRNA genes and regulatory RNA structure on mRNAs. We discuss the strengths and weaknesses of the different strategies and how they can complement each other

Mapping the RNA Chaperone Activity of the T. brucei Editosome Using SHAPE Chemical Probing

Mitochondrial pre-mRNAs in African trypanosomes adopt intricately folded, highly stable 2D and 3D structures. The RNA molecules are substrates of a U-nucleotide-specific insertion/deletion-type RNA editing reaction, which is catalyzed by a 0.8 MDa protein complex known as the editosome. RNA binding to the editosome is followed by a chaperone-mediated RNA remodeling reaction. The reaction increases the dynamic of specifically U-nucleotides to lower their base-pairing probability and as a consequence generates a simplified RNA folding landscape that is critical for the progression of the editing reaction cycle. Here we describe a chemical mapping method to quantitatively monitor the chaperone-driven structural changes of pre-edited mRNAs upon editosome binding. The method is known as selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE). SHAPE is based on the differential electrophilic modification of ribose 2'-hydroxyl groups in structurally constraint (double-stranded) versus structurally unconstrained (single-stranded) nucleotides. Electrophilic anhydrides such as 1-methyl-7-nitroisatoic anhydride are used as probing reagents, and the ribose 2'-modified nucleotides are mapped as abortive cDNA synthesis products. As a result, SHAPE allows the identification of all single-stranded and base-paired regions in a given RNA, and the data are used to compute experimentally derived RNA 2D structures. A side-by-side comparison of the RNA 2D folds in the pre- and post-chaperone states finally maps the chaperone-induced dynamic of the different pre-mRNAs with single-nucleotide resolution

Genome-wide mapping of RNA structure using nuclease digestion and high-throughput sequencing.

RNA structure is important for RNA function and regulation, and there is growing interest in determining the RNA structure of many transcripts. Here we provide a detailed protocol for the parallel analysis of RNA structure (PARS) for probing RNA secondary structures genome-wide. In this method, enzymatic footprinting is coupled to high-throughput sequencing to provide secondary structure data for thousands of RNAs simultaneously. The entire experimental protocol takes ∼5 d to complete, and sequencing and data analysis take an additional 6-8 d. PARS was developed using the yeast genome as proof of principle, but its approach should be applicable to probing RNA structures from different transcriptomes and structural dynamics under diverse solution conditions. Nat Protoc 2013 May; 8(5):849-69

Continuously tunable nucleic acid hybridization probes

In silico–designed nucleic acid probes and primers often do not achieve favorable specificity and sensitivity tradeoffs on the first try, and iterative empirical sequence-based optimization is needed, particularly in multiplexed assays. We present a novel, on-the-fly method of tuning probe affinity and selectivity by adjusting the stoichiometry of auxiliary species, which allows for independent and decoupled adjustment of the hybridization yield for different probes in multiplexed assays. Using this method, we achieved near-continuous tuning of probe effective free energy. To demonstrate our approach, we enforced uniform capture efficiency of 31 DNA molecules (GC content, 0–100%), maximized the signal difference for 11 pairs of single-nucleotide variants and performed tunable hybrid capture of mRNA from total RNA. Using the Nanostring nCounter platform, we applied stoichiometric tuning to simultaneously adjust yields for a 24-plex assay, and we show multiplexed quantitation of RNA sequences and variants from formalin-fixed, paraffin-embedded samples